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Geodesic motion determines important features of spacetimes. Null unstable geodesies are closely 
related to the appearance of compact objects to external observers and have been associated with 
the characteristic modes of black holes. By computing the Lyapunov exponent, which is the inverse 
of the instability timescale associated with this geodesic motion, we show that, in the eikonal 
limit, quasinormal modes of black holes in any dimensions are determined by the parameters of 
the circular null geodesies. This result is independent of the field equations and only assumes a 
stationary, spherically symmetric and asymptotically flat line element, but it does not seem to be 
easily extendable to anti-de Sitter spacetimes. We further show that (i) in spacetime dimensions 
greater than four, equatorial circular timelike geodesies in a Myers-Perry black hole background are 
unstable, and (ii) the instability timescale of equatorial null geodesies in Myers-Perry spacetimes 
has a local minimum for spacetimes of dimension d > 6. 

PACS numbers: 04.70.Bw,04.50.Gh,05.45.-a 



I. INTRODUCTION 

Geodesies in black hole spacetimes have been ex- 
tensively studied, both in four and higher dimensional 
spacetimes, with and without a cosmological constant. 
Geodesies may display a rich structure and they con- 
vey important information on the background geometry. 
Among the different kinds of geodesic motion, circular 
geodesies are especially interesting. For instance, the 
binding energy of the last stable circular timelike geodesic 
in the Kerr geometry is related to the gravitational bind- 
ing energy that can be radiated to infinity, and it can 
be used to estimate the spin of astrophysical black holes 
through observations of accretion disks [H, & 9 • 

It was shown many years ago that null geodesies also 
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play an important role. The optical appearance of a star 
undergoing gravitational collapse depends crucially on 
the circular unstable null geodesic, which also explains 
an exponential fade-out of the collapsing star's luminos- 
ity Null geodesies are also very useful to explain 
the characteristic modes of a black hole - the so-called 
quasinormal modes (QNMs) [1, 0|- These "free" modes 
of vibration can be interpreted in terms of null particles 
trapped at the unstable circular orbit and slowly leaking 
out [1, SIM [II El- The real part of the complex QNM 
frequencies is determined by the angular velocity at the 
unstable null geodesic; the imaginary part is related to 
the instability timescale of the orbit, a quantity which is 
seldom considered in geodesic studies, with some note- 
worthy exceptions (see e.g. QJ, [lj, GEES d El EH)- 
Furthermore, there is some evidence [l5[ that unstable 
circular orbits could yield information on phenomena oc- 
curring at the threshold of black hole formation in the 
high-energy scattering of black holes, a process of interest 
in fundamental physics for a variety of reasons [H HH • 

In this work we clarify some aspects of the relation be- 
tween unstable null geodesies, Lyapunov exponents and 



2 



quasinormal modes. In Section [II] we derive a simple for- 
mula for the principal Lyapunov exponent A in terms of 
the second derivative of the effective potential for radial 
motion V r : 



A 



(1) 



where t is coordinate time. Throughout this work, a dot 
denotes a derivative with respect to proper time and a 
prime stands for derivative with respect to areal radius 
r. The result above is valid for a wide class of spacetimes 
and geodesies, including stationary spherically symmet- 
ric spacetimes and equatorial orbits in the geometry of 
higher-dimensional, rotating (Myers-Perry) black hole so- 
lutions. 

In Sect ion ITU we show that the relation between QNMs 
and unstable circular null geodesies is quite general, be- 
ing valid in the eikonal limit for any static, spherically 
symmetric, asymptotically flat spacetime. More specifi- 
cally, we show that the angular velocity £l c at the unsta- 
ble null geodesic and the Lyapunov exponent, determin- 
ing the instability timescale of the orbit (see for instance 
[la , IT4T | ) agree with analytic WKB approximations for 
QNMs IIH: 



wqnm =il c l — i(n+ 1/2) |A| 



(2) 



where n is the overtone number and / is the angular 
momentum of the perturbation. The WKB results are 
formally valid only in the eikonal regime (I 1), but 
they seem to yield surprisingly accurate predictions even 
for low values of I fl2l [25j . A simple derivation of the 
Lyapunov exponent for spherically symmetric, asymptot- 
ically flat spacetimes, patterned after the original QNM 
calculation by Ferrari and Mashhoon [l^, EH, is pro- 
vided in Appendix [XJ For the important case of a d- 
dimensional Schwarzschild- Tangherlini (2(| black hole so- 
lution we find that the critical exponent defined by Pre- 
forms and Khurana [l5| can be determined analytically 
to be 
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By exploring the connection between QNMs and null 
geodesies, we also find a simple analytical result for the 
quasinormal frequencies of near-extremal Schwarzschild- 
de Sitter black holes in d = 4: 



wqnm = k+ [I - i (n + 1/2)] 



(4) 



where k + denotes the surface gravity. In the eikonal 
limit, the above result agrees with that found in [271 ]. 

In Section lTVl we analyze the higher-dimensional rotat- 
ing black hole solutions found by Myers and Perry [2^ |. 
In d = 5 we can compute A analytically. The Lyapunov 
exponent goes to zero as one approaches extremality in 
d — 4, 5 spacetime dimensions. However, no such be- 
havior is observed for d > 5: the Lyapunov exponent 



(normalized by the orbital frequency) has a local mini- 
mum, which may be related to a possible instability of 
the system first suggested by Emparan and Myers [23 ]. 
In Appendix [B] wc study in some detail timelike circular 
geodesies in the Myers-Perry spacetime and show that 
equatorial circular orbits are always unstable for d > 4. 
Finally, in Appendix [C] we discuss issues in generalizing 
our results to non- asymptotically flat spacetimes. 



II. LYAPUNOV EXPONENTS AND GEODESIC 
STABILITY 

Lyapunov exponents are a measure of the average rate 
at which nearby trajectories converge or diverge in the 
phase space. A positive Lyapunov exponent indicates a 
divergence between nearby trajectories, i.e., a high sensi- 
tivity to initial conditions. A geodesic stability analysis 
in terms of Lyapunov exponents begins with the equa- 
tions of motion schematically written as 



dXj 
dt 



and linearized about a certain orbit: 
dSXi(t) 



dt 



Here 



K i:i (t) = 



dH, 



(5) 



(G) 



(7) 



x,(t) 



is the linear stability matrix E3- The solution to the 
linearized equations can be written as 

SK i (t)=L i j(t)5X j (0) (8) 

in terms of the evolution matrix L^-(t), which must obey 



(t) 



(9) 



and Lij(fi) 



,j . A determination of the eigenvalues of 
Lij leads to the principal Lyapunov exponent A, which is 
the quantity we focus on. Specifically 



A = lim — log ' ^ 



t 



(0 



(10) 



We now restrict attention to a class of problems for 
which one has a two dimensional phase space of the form 
Xi(t) = (p r ,r). This includes circular orbits in station- 
ary spherically symmetric spacetimes and equatorial cir- 
cular orbits in stationary spacetimes, such as the Myers- 
Perry metric considered in Section HVl Linearizing the 
equations of motion with Xi(t) = (j> r ,r) about orbits of 
constant r we get 
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(11) 
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where 



d 
dr 



(ig rl 



X SC 



(12) 
(13) 



and C is the Lagrangian for geodesic motion (see below 
for explicit examples). Therefore, for circular orbits, the 
principal Lyapunov exponents can be expressed as 



X = ±^/K 1 K 2 . (14) 

From now on we will drop the ± sign, and simply refer to 
the "Lyapunov exponent" . From the equations of motion 
it follows that 

d_S£ _ SC 
dr Sr Sr ' 

and 

d 



(15) 



d SC d 
dr Sr dr 
Using the definition of V r , 



dr 



= V r 



(16) 



(17) 



we can rewrite this as 
SC _ 
~Sr~ ~ ' 



2g r 

For circular geodesies V r - 
reduces to 

A = 



dr 
V 



{glrVr) ■ 



o tec 



(18) 

and Eq. (TJD 




(19) 

we can define a 



(20) 



Following Pretorius and Khurana [If 
critical exponent 

n c _ t x 

7 ~ 2ttA T n ' 

where we introduced a typical orbital timescale Tq = 
2tt/H, c and an instability timescale T\ = 1/A (note that 
in Ref. the authors use a different definition of the 
orbital timescale, Tq = 2ir/p, with ip an angular coordi- 
nate). Then we get 



1 
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III. SPHERICALLY SYMMETRIC 
SPACETIMES 

We will consider a stationary spherically symmetric 
background 

ds 2 = f{r)dt 2 - -^-dr 2 - r 2 dfl 2 d _ 2 , (22) 
9{r) 

where f(r) and g(r) are functions to be determined by 
solving the field equations, dil 2 l _ 2 is the metric of the 
(d - 2)-sphere and A d _ 2 = 2ir (d - 1 ^ 2 /T[(d - l)/2] is the 
area of the unit (d — 2)-sphere. Since we do not spec- 
ify the field equations, our results are valid for any field 
equations admitting spherically symmetric, asymptoti- 
cally flat solutions. The last property will be required 
to enforce the correct boundary conditions in the WKB 
calculations of Section IIII CI 



A. Circular orbits 

To compute the geodesies in the spacetime (|22|) we 
follow Chandrasekhar (34|. Let us restrict attention to 
equatorial orbits, for which the Lagrangian is 



2C = f(r) i 2 



1 



9(r) 
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r — r ip . 



(23) 



where ip is an angular coordinate. The generalized mo- 
menta derived from this Lagrangian are 



p t = f(r) i=E= const , 
p v = — r~ ip = —L = const , 



Vr = 



-r 2 ip \ 
1 

'W) 



(24) 
(25) 

(26) 



The Lagrangian is independent of both t and ip, so it 
follows immediately that p t and p v are two integrals of 
motion. Solving ([24|) -([25 |l for t,<fi we get 



E 

W) 



(27) 



The Hamiltonian is given by 



For circular null geodesies in many spacctimcs of inter- 
est V" > 0, which implies instability. A quantitative 
characterization of this instability can be achieved by 
computing the timescale associated with it. In most of 
this pap er we will use the method outlined above (see 
also 113. l3ll. |32| 1. but there are alternative approaches 
EElSMEIa]. In Appendix O for example, we recover 
the results of the next Section following a stability anal- 
ysis due to Ferrari and Mashhoon [l(| [HI • 

The discussion above is still rather general, assuming 
only that the variables in the equations of motion form a 
two-dimensional plane in phase space. We now specialize 
to spherically symmetric spacetimes. 



2TL = 2 yp t t + p v p + p r f — C 
1 .0 



f(r)t 2 



g( r ) 



r 2 <p 2 



1 



Et — Lip — — -f = 81 = const . 

9\r) 



(28) 



Here S\ — 1,0 for time-like and null geodesies, respec- 
tively. Inserting Eq. 127)1 m (EE]) an d using the definition 
(QUI we get 



V r = g(r) 



E 2 
f(r) 



Si 



(29) 
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1. Timelike geodesies 
The requirement V r = = for circular orbits yields 

(30) 



2 2/ 2 2 r 3 /' 

£ 2 = L — J 



2/ - rf 



where here and below all quantities are evaluated at the 
radius of a circular timelike orbit. Since the energy must 
be real, we require 

2/-r/'>0. (31) 
The second derivative of the potential is 

and the orbital angular velocity is given by 

n = i = (L\'\ (33) 



t \2r 



2. Null geodesies 



Bearing in mind that 2f — rf > and that unstable 
orbits are defined by V" > 0, we can see that A will be 
real whenever the orbit is unstable, as expected. In d = 4 
this formula reduces to well-known results [3]. We also 
get 



2vrA 



= (2-) 2 



-3ff + 2r^/'-r^ r 



(39) 



2. Null geodesies 



Using Eqs. (HU), ([27]), (O and (JM]), for circular null 
geodesies we find 



* / T—T. 



(40) 

In the last equality we made use of ([33]) and we defined 
the "tortoise" coordinate 



(41) 



Circular null geodesies satisfy the conditions: 



E 

T 



2/c 




(34) 
(35) 



Here and below a subscript c means that the quantity in 
question is evaluated at the radius r — r c of a circular null 
geodesic. An inspection of (|3"5"|) shows that circular null 
geodesies can be seen as the innermost circular timelike 
geodesies. In this case 



K'(r c ) 



L 2 9c 
rife 



and the coordinate angular velocity is 



n c = % = 
t 



A 

2r r 



1/2 



c l/2 



(36) 



(37) 



B. Lyapunov exponents 

1. Timelike geodesies 

Using Eqs. (19]), (27]) and (J32J) to evaluate the Lya- 
punov exponent at the circular timelike geodesies, we get 



(38) 




y(2/-r/')K"W 



C. Unstable null geodesies and quasinormal modes: 
comparison with WKB results 

WKB methods [H \M provide an accurate 
approximation of QNM frequencies in the eikonal limit 
for spacetimes where the wave equation can be cast in 
the form 



^-2* + Qo* = 0, 



(42) 



where r» is a convenient "tortoise" coordinate, ranging 
from — oo to +oo. In particular, one gets the QNM con- 
dition 



)(ro) 



t(n + l/2). 



(43) 



where = d 2 Q /drl and Eq. (|4*3"]) is evaluated at the 
extremum of Qq, i.e. the point tq at which dQo/dr* = 0. 
We note that this result is strictly valid only for asymp- 
totically flat spacetimes, or for spacetimes admitting 
wavelike solutions at spatial infinity. It is not valid for 
anti-de Sitter (AdS) backgrounds. In a spacetime of the 
form (|2"2"j) . we find that the Klein-Gordon equation can be 
written as in Eq. (|4"2"j) with the tortoise coordinate (|41[) . 
In the eikonal limit (Z — > oo) we get 



(44) 



It is known that scalar, electromagnetic and gravitational 
perturbations of static black holes in higher dimensions 
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have the same behavior in the eikonal limit [35J, [36J, [37j . 
In other words, there is a well-defined geometric-optics 
(eikonal) limit where the potential for a wide class of 
massless perturbations is "universal" . For Qo above we 
find that the extremum of Qo satisfies 2/(ro) = fof'^fo), 
i.e. ro coincides with the location of the null circular 
geodesic ro = r c , as given by Eq. l|35p . Furthermore, the 
WKB formula (|4"3"| allows one to conclude that, in the 
large-/ limit, 



WQNM 




fc 



<p_ 

dr 2 



f 



(45) 



Comparing with Eqs. ((37)) and ([40]) we find that 

wqnm =n c l- »(n+ 1/2) |A|. (46) 

This is one of the main results of this paper: in the 
eikonal approximation, the real and imaginary parts of 
the QNMs of any spherically symmetric, asymptotically 
flat spacetime are given by (multiples of) the frequency 
and instability timescale of the unstable circular null 
geodesies. 



The calculation of the Pretorius-Khurana [if 
ponent yields 
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2ttX 



Tx 



1 



2n^d - 3 



critical ex- 



(50) 



where in the last equality we made use of Eq. (|49|) . This 
result is in excellent agreement with numerical calcula- 
tions by Merrick and Pretorius [39]. Small values of 7 
correspond to a strong Lyapunov instability, so the insta- 
bility is more pronounced for large spacetime dimensions. 
By relating the geodesic orbital frequency and instability 
timescale to the QNM frequencies in the eikonal limit, 
we get 



WQNM 



Z-n/d-3(ra + l/2). 



(51) 



This is in complete agreement with known analytical and 
numerical results in d — 4 (40| and higher dimensions 
HHEl, 43], for asymptotically flat spacetimes. 



Near-extremal Schwarzs child- de Sitter spacetime in four 
dimensions 



1. Higher- dimensional Schwarzschild black hole 

Let us consider a more specific example: the higher- 
dimensional Schwarzschild- Tangherlini metric [26j 



ds 2 



fd^-f^dr 2 



l dn 



d-2 1 



f(r) 



1 



d-3 



which includes the well-known four-dimensional geom- 
etry as a special case. Here dfl 2 _ 2 is the metric of 
the (d — 2)-sphere, and the horizon radius r + is re- 
lated to the mass M of the spacetime via M = (d — 
2) Ad- 2?+ 3 /(16tt). For timelike geodesies we find that 
the orbits must satisfy 



r > r c 



(48) 



where r c is the radius of the circular null geodesic. With 
the requirement (f3"Tj) we have V" > for all d > 4, 
and therefore all circular orbits are unstable for d > 4 
[26L [38| . The four-dimensional case is special: one gets 
F r " = 4Mr(6M - r)/(2r - 6M). Thus, in four spacetime 
dimensions there are stable circular orbits for any r > 
6M. The circular orbits with radius 3M < r < 6M are 
all unstable. For light-like geodesies, one has V"(r c ) = 
L 2 {2d — 6) r^ 4 > 0. Therefore circular null geodesies are 
always unstable for d > 4. The angular velocity at r c is 
given by 



d-3 



1 



(49) 



A non-trivial example concerns a non-asymptotically 
flat spacetime, the near-extremal Schwarzschild-de Sit- 
ter (SdS) spacetime in four dimensions. General SdS 
spacetimes have a metric of the form (p?!?|) with f(r) — 
g{r) = 1 - 2M/r - (r/L ds ) 2 . M denotes the black-hole 
mass and is the de Sitter curvature radius, related to 
the cosmological constant A by L^ s — 3/ A. The space- 
time possesses two horizons: the black-hole horizon is at 
r = r + and the cosmological horizon is at r = rc , where 
r Co > r + . The function / has zeroes at r+, rc , and 
ro = — (r+ + rco)- In terms of these quantities, / can be 
expressed as 



/ 



1 



(r - r+)(r Co - r)(r - r ). 



(52) 



It is useful to regard r + and rco as the two fundamental 
parameters of the SdS spacetime, and to express M and 



L 2 S as functions of these variables: 



1 1 — 



2 + r + r Co + r Co 2 ■ 



2ML 2 ds = r + r Co (r+ + r Co ) . 



(53) 
(54) 



We also introduce the surface gravity k+ associated with 



the black hole horizon r 
plicitly, we have 



r + : k + = \{df/dr) r -_ 



K + 



(r Co - r + )(r + - r ) 



Ex- 



(55) 



Let us now specialize to the near-extremal SdS black 
hole, which is defined as the spacetime for which the cos- 
mological horizon rco is very close (in the r coordinate) 
to the black hole horizon r + , i.e. 



rco - r + < r 4 



(56) 



6 



For this spacetime one can make the following approxi- 
mations: 



where 



r ~ -2r\ ■ L 2 S ~ 3r 2 ; M ~ -± ; k h 



2r 2 



(57) 

Note that k+ is infinitcsimally small, since we assume 
Eq. ([56]) holds. For null geodesies, an exact solution can 
be found with 



(58) 



The angular velocity at this radius is given by 



^s^ = -l + i(l-^) _ f . (59) 
and reduces to 



r 2 + a 2 cos 2 i? . 



i = r 2 + a 2 - //r 5 ~ d , (64) 

and dCli^ denotes the standard metric of the unit (d — 
4)-sphere [28]. This metric describes a rotating black 
hole in an asymptotically flat, vacuum spacetime with 
mass and angular momentum proportional to \i and \xa, 
respectively. Hereafter we assume fi > and a > 0. 

The event horizon is located at r = r + such that 
A| r=r+ =0. In the standard four-dimensional case, an 
event horizon exists only for a < ji/2. In d = 5 an event 
horizon exists only when a < J]i, and it shrinks to zero 
area in the extreme limit a — > JJi. On the other hand, 
when d > 6, which is the part of the parameter space 
that we shall focus on, A = has exactly one positive 
root for arbitrary a > 0. This means there is no bound 
on a, or (loosely speaking) that there are no extremal 
Kerr black holes in higher dimensions. 



_ (r Co - r+) _ 
Uc ~ 2r 2 ~ 4 



(60) 



A. Circular geodesies in the equatorial plane 



in the near-extremal regime. As in the case of four- 
dimensional Schwarzschild black holes, Eq. ([50)) , we find 



A/fi c =l. (61) 

This formula predicts 

wqnm = k+ [I - i (n + 1/2)] . (62) 

The QNMs of this spacetime are known in closed form 
[13] and agree with (|62|) in the eikonal limit. 



IV. MYERS-PERRY BLACK HOLES 

In four dimensions there is only one possible rotation 
axis for an axisymmetric spacetime, and there is there- 
fore only one angular momentum parameter. In higher 
dimensions there are several choices of rotation axis and 
there is a multitude of angular momentum parameters, 
each referring to a particular rotation plane. Rotating 
black hole solutions in higher dimensions are known as 
Myers-Perry black holes [28|]. We focus on the simplest 
case for which there is only one angular momentum pa- 
rameter, that we shall denote by a. The metric of a 
d-dimensional Myers-Perry black hole with only one non- 
zero angular momentum parameter in Boyer-Lindquist- 
type coordinates is given by (here we adopt the notation 
commonly used in related works [H, H EE HI, H3 ) , 

l2 A - a 2 sin 2 , 2 2a(r 2 + a 2 - A) sin 2 , , 
ds 2 = dt 2 + — i '- dtdip 

(r 2 + a 2 ) 2 - Aa 2 sin 2 # . 



A 



sin §dip 

Zj 

dr 2 - Erfi? 2 - r 2 cos 2 MVL 2 ^, 



To write down the geodesic equations in the Myers- 
Perry spacetime we follow Chandrasekhar (34|. Let us 
restrict attention to orbits in the equatorial plane (i? = 
0,i? = it/2), for which the appropriate Lagrangian is 

2C = g tt i 2 + 2g tv itp + g rr r 2 + g w Cp 2 . (65) 

The generalized momenta following from this Lagrangian 
are 



Pi 

Pr 



gut + gt v <P = 

grrf • 



E = const , 
= — L = const . 



(66) 
(67) 
(68) 



The Lagrangian is independent of both t and ip, so it 
follows immediately that p t and p v are two integrals of 
motion. Solving (f6"6"|) - (f6"T)) for i , (p we find 



1 

v = A 
1 

1 ^ A 



E 



1 - 



,d-3 



L 



a 2 fi 

r d-3 



£7- 



The Hamiltonian is given by 

2TL = 2 (p t i + p v <f >+ p r f — C 



d-3 



L 



(69) 
(70) 



= 1- 



■*d— 3 



■ 2aa . 



A 



r 2 



o 2 /z 

r d— 3 



Et- Lip 



const . 



(71) 



Here <5i = 1 , for time-like and null geodesies, respec- 
tively. Inserting Eqs. (|6l? ]) -([70" ]) in Eq. ([7TJ) we get 

(72) 



(63) r 2 V r 



V r , 

r 2 E 2 



r d—3 



(aE-L) 2 + (a 2 E 2 -L 2 )-6 1 A 
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In d = 4, recalling that /i = 2M, we recover the well- 
known results for the Kerr geometry [34| . In five dimen- 
sions we recover the results by Frolov and Stojkovic 
if we specialize their equations to only one rotation pa- 
rameter. 

The conditions for the existence of circular geodesies, 
V r = V{ = 0, lead to the following equations: 

= r 2 E 2 +^r a - d (aE-L) 2 + (a 2 E 2 -L 2 )-S 1 A, 
= Ar 2 E 2 — (d — 5)[ir 3 ~ d (aE — L) 2 + 2(a 2 E 2 — L 2 ) 
-5 1 (2A + rA'). (73) 

Eliminating the term a 2 E 2 — L 2 we get 



2rE* 



(d - 3)Lir 2 - d (aE - L) 2 - A'<*i = 0. 



(74) 



result [34[ to the general case of Myers-Perry spacetimes. 
It would be interesting to investigate whether or not this 
is a general property of any stationary spacetime. 

For corotating orbits with a ^/(^-s) and d > 5, the 
following analytical approximations are valid: 



D r — a 




A*(d-3) \2aV2 V d-3, 
(corotating, a — > oo) . (82) 



In particular, when the rotation is very large the radius 
of the corotating orbit is at fixed relative distance from 
the horizon: 



B. Circular null geodesies 



For light-like geodesies (Si 
conditions 



0) we get the explicit 



r 2 c E 2 + (ir 3 c - d (aE~L) 2 + (a 2 E 2 
2r c E 2 - (d - 3)r^ _d n(aE - L) 2 . 



L 2 ), (75) 
(76) 



The above equations can be simplified by the introduc- 
tion of the impact parameter D c = L/E: 



= r 2 c +[ir 3 c - d (a-D c ) 2 + (a 2 -D 2 ) 

.2— d/„ n \2 



= 2r c -{d-3) fJ ,r 2 c - a (a-D c ) 
From Eq. (JTSJ) we get 



2r, 



d-l 



(d-3)ti 



(77) 
(78) 



(79) 



Notice that Eq. ([77]) is satisfied if and only if \D C \ > a. 
For counterrotating orbits, we have \D C — a\ = —{D c — a) 
and this case corresponds to the upper sign in the above 
equation, while for corotating orbits, \D C — a\ = +(D c —a) 
and this case corresponds to the lower sign in Eq. (fT9"|) . 

Substituting Eq. ([79")) in ([77]) . we find an equation for 
the radius of circular null geodesies: 



d-l 
d-3 



ri ± 2a\ 



2r, 



d-l 



2r 



d-l 



(d — 3)/i (d — 3)/i 



= 0. 



(80) 



In d = 4 we recover the well-known results [30j . 

An important quantity for the analysis of the null 
geodesies is the angular frequency at the null geodesic 
0,: 



r+ 



d-l 



2^2(d-3) 



(83) 



This may be explained by the fact that the angular ve- 



locity of the horizon, Q r 



rj_-{-a- 



■ , also has a maximum 



and then decreases to zero for very large a. 

On the other hand, counterrotating orbits are well de- 
scribed by 



r (d-i)/2 M ay / 2f i(d-3), 
1 

(counterrotating , a — * oo) . (84) 



We can express the radius of counterrotating orbits in 
terms of the horizon radius in the limit of very large 
rotation: 



^^/2(rf^3) 



(85) 



Notice how counterrotating orbits must be located very 
far away from the horizon as rotation increases. This is a 
consequence of having a "strong" ergoregion, extending 
through a large region in space. 

For d = 5 the previous equations simplify considerably, 
and we can find a simple solution: 



r c = \/2yV±a^7i, 

— = n f= 1 

D c -a=F2yf]I 



(86) 
(87) 



3-d 



+ (1 - ^r 3 c - d )D c 



{rl + a 2 + a Ve ) 



afirl d D c 



h (81) 



where we have used Eqs. (|T9"|) and (T5U|) . Therefore the fre- 
quency of equatorial null geodesies is the inverse of their 
impact parameter. This generalizes the four-dimensional 



C. Circular timelike geodesies 

Timclike geodesies in a Myers-Perry spacetime are 
studied in Appendix [5] In Section IB II we show that 
for d > 4, there are no stable equatorial circular orbits 
in this spacetime. This extends the instability proof by 
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FIG. 1: Ratio of the timelike orbital frequency fi to the orbital frequency fl c of corotating (left) and counterrotating (right) 
null geodesies as a function of r/r c , with r c the radius of the null geodesic. 



Tangherlini [26| to rotating black holes, and the instabil- 
ity proof by Frolov and Stojkovic [48] to a general number 
of spacetime dimensions. 

The energy and angular momentum of timelike circu- 
lar geodesies are studied in Section IB 2i where we also 
consider the orbital frequency f2 of general circular time- 
like geodesies. In fact, it is possible to obtain a simple 
expression for the ratio f2/f2 c [see Eq. (IB25[) in Appendix 

EI- 



n 
n7 



This quantity is plotted in Fig. Q] for /i = 2 and selected 
values of a and d. There is clearly a change in behavior 
for d > 5 and large rotation, exemplified here for d = 6. 
This will be explored in more detail in the next Section. 

D. Lyapunov exponents 

Applying Eq. (|19p to the case of Myers-Perry black 
holes we find 



(89) 



Using Eq. (|86|) we can analytically compute the d = 5 
case, and we get 



A c 



intcr ^7 = v/2 
A co V/I = V2 



a — 2y / /I 
a + 



(90) 
(91) 



The general case requires a numerical treatment. In 
Fig. fj] we show the Lyapunov exponent normalized by 
the angular velocity, A/f2 c , as a function of rotation a. 
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FIG. 2: Dimensionless instability exponents A/f2 c as a func- 
tion of rotation for several spacetime dimensions d. We use 
units such that fi = 2. Solid lines refer to corotating orbits, 
dashed lines to counterrotating orbits. 



For corotating geodesies in d = 4 and d = 5, A asymp- 
totes to zero near extremality (a — > /x/2 and a — ► ^//I, 
respectively). This does not happen for d > 5 and large 
rotation parameters: in this case, A/f2 c has a local min- 
imum. This is consistent with QNM calculations in five 
[4^| and higher dimensions [43| . The local minimum may 
be related to a change in behavior correspondin g t o a 
black hole — > black brane transition, observed in [29| in 
relation with a conjectured instability of these systems 
for very large rotation rates. 

The locations of minima in A/f2 c are given in Table |T] 
In this Table we also list a quantity considered in 29] as 
a possible indicator of a "transition" point, in this case 
the rotation at which the temperature has a minimum. 
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TABLE I: The transition point atrans- The second column 
refers to the critical value of rotation a for which the corotat- 
ing Lyapunov exponent has a minimum. The third column 
refers to the point at which the temperature has a minimum, 
an indicator considered in [29l ]. 

d Lyapunov Temperature 



6 


1.15 


1.37 


7 


0.84 


1.28 


8 


0.64 


1.22 


50 


0.05 


1.02 



The two quantities are roughly consistent for small d, but 
not for large d. 



V. CONCLUSIONS AND FUTURE WORK 

We have shown that for all spherically symmet- 
ric spacetimes, in a geometrical optics approximation, 
QNMs can be interpreted as particles trapped at unsta- 
ble circular null geodesies and slowly leaking out. The 
leaking timescale is given by the principal Lyapunov ex- 
ponent, for which we obtained a fairly simple expression, 
Eq. (HHJ), in terms of the second derivative of the effective 
radial potential for geodesic motion. This simple, intu- 
itive relation between QNMs and circular null geodesies 
is valid for all asymptotically flat, spherically symmetric 
black hole spacetimes. 

Some aspects of our investigation deserve further anal- 
ysis. The interpretation of QNMs in terms of unsta- 
ble circular null geodesies is valid in all generality only 
for spherically symmetric, asymptotically flat spacetimes. 
Once we break the azimuthal degeneracy, things get a bit 
more complex. For instance, in d — 4 it is known that 
equatorial geodesies can account for the I = \m\ modes 
of Kerr-Newman black holes [HI HH, l5(ij |. but it is un- 
clear whether this analogy can be extended to modes 
with I ^ \m\. Perhaps modes with I m can be ex- 
plained in terms of more general (e.g., non-equatorial) 
geodesies. Besides improving our intuitive physical un- 
derstanding of ringdown radiation, a deeper exploration 
of this analogy could have important implications for the 
interpretation of numerical simulations of black hole bi- 
nary mergers and their use in gravitational-wave data 
analysis [Bll. I52L f53j] . 

Another limitation of our results concerns their exten- 
sion to non asymptotically flat (e.g., AdS) backgrounds. 
In the large damping limit, a certain class ofQNMs has 
been associated with radial geodesies [54| Ha HfJ- It 
would be very interesting to extend this analysis to the 
eikonal (large-/) limit. A possible starting point could be 
the 2 + 1 dimensional Bahados-Teitelboim-Zanelli (BTZ) 
black hole 15711 . for which QNM frequencies are known 
analytically [581 ]. Quite apart from the geodesic anal- 
ogy, the large-/ limit is interesting per se. It turns out 
that the imaginary part of Schwarzschild-anti-de Sitter 



(SAdS) QNMs decreases with / [56|. Thus, if excited con- 
siderably, large-/ modes could dominate the black hole's 
response to perturbations. A more thorough investiga- 
tion of the large-/ limit of QNMs is necessary. 

Interesting physical phenomena could occur in (hypo- 
thetical) spacetimes for which timelike circular geodesies 
have a frequency equal to (or larger than) the frequency 
of unstable null geodesies. This would raise the inter- 
esting possibility of exciting QNMs by orbiting parti- 
cles, possibly leading to instabilities of the spacetime. 
It would be interesting to find general conditions un- 
der which spacetimes possess stable null geodesies; stable 
null geodesies may also be associated with instability (or 
marginal stability) of spacetimes. 
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APPENDIX A: THE INSTABILITY TIMESCALE 
OF CIRCULAR GEODESICS: A SIMPLER 
DERIVATION 

Perhaps the simplest way to determine the instabil- 
ity timescale associated with circular null geodesies is 
through a consideration of the equations defining these 
geodesies d, [13, EU- Indeed, this was the approach 
originally adopted by Ferrari and Mashhoon to compute 
quasinormal modes of Kerr-Newman black holes in the 
eikonal limit [lol . [TlT | . In this Appendix we will rederive 
the instability parameter A within their approach. 

Consider small perturbations of a bundle of test null 
rays in the unstable equatorial circular orbit around a 
black hole described by the metric (f2"2")l . First, rescale 
the affine parameter s so as to be the coordinate time 
t, and consider the following values of the unperturbed 
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geodesic: 



t = s , r = 77. 



7T 

2 



= 0= 



(Al) 



The slightly perturbed equatorial null orbit is given by 
s = t + eb(t), (A2) 
r = r c (l + eh(t)), (A3) 
= tt/2, (A4) 
= n c (t + efe(t)), (A5) 

where |e| <C 1 denotes the dimensionless amplitude of 
the perturbation. Considering the leading terms in Eqs. 
(f2"T)) - (|29[) . together with the boundary conditions that 
b(t), h(t) and k(t) vanish at t = 0, yields 



h(t) 


= sinh(At) , 


(A6) 


b(t) 


= ^(cosh(At)-l), 


(A7) 


k(t) 


= 0. 


(A8) 



The parameter A characterizes the decay rate and is de- 
termined by 

X " = j2 y >c) = ^KVc), (A9) 

This result agrees with Eq. (j4"0)) , obtained through the 
principal Lyapunov exponent. 

Finally, the instability timescale can also be derived (at 
least in the usual Schwarzschild geometry) by considering 
a special class of geodesies: inspiralling geodesies that 
asymptote to the light ring when t — > 00. These geodesies 
are considered in Chandrasekhar's book [34], and it is 
straightforward to compute how they approach the light 
ring. 



APPENDIX B: TIMELIKE GEODESICS IN THE 
EQUATORIAL PLANE OF MYERS-PERRY 
SPACETIMES 

In a d-dimensional Myers-Perry black hole with only 
one non-zero angular momentum parameter, the radial 
equation for geodesies in the equatorial plane can be cast 
in the form r 4 f 2 = V with V = r 4 V r as given in Eq. (fT2")l . 
The first derivative of V with respect to r is given by 

V = 4r 3 E 2 - (d-5)fir 4 - d (aE - L) 2 + 2r(a 2 E 2 - L 2 ) 

-Si [4r 3 + 2ra 2 + (d - 7)fir 6 ~ d ] . (Bl) 

The conditions for the existence of circular orbits are 

V = and V = 0: 







r 2 E 2 + f ir 3 - d (aE-L) 2 
+ (a 2 E 2 - L 2 )-(5iA, (B2) 
= Ar 2 E 2 — (d — 5)fir 3 ~ d (aE - L) 2 + 2(a 2 E 2 - L 2 ) 
- 5 l (2A + rA') . (B3) 

Eliminating the term a 2 E 2 — L 2 one finds 



2rE z 



(d — 3)[ir 2 ~ (aE - L) 2 - A' Si = . (B4) 



1. The instability of geodesies in the equatorial 
plane of Myers-Perry spacetimes 



The equations (|B2j) and (IB4j) will be used here to ob- 
tain the values of energy E and angular momentum L 
associated to circular timelike (Si = 1) geodesies. Using 
(|B4j) and introducing the new quantities x — L — aE, 
M = /Lt/2, and the reciprocal radius u = l/r, we obtain 
an expression for E 2 : 



E 2 = [1 - (5 - d)Mu d - 3 ] +(d- 3)Mu d 
With this expression for E 2 , Eq. (|B2[) leads to 
2aExu = x 2 [(d- l)Mu d - 3 ~ l] u 



1 x 2 . 



[a 2 u -(d- 3)Mu 



d-41 



(B5) 



(B6) 



For d — 4 spacetime dimensions, the above equations are 
identical to those obtained by Chandrasekhar [3J|. We 
can now eliminate E in equations (|B5[) and (|B6[) to obtain 
a quadratic equation for x 2 : 



= x A u 2 



-2x 2 u 



[(d - l)Mu d -' 3 - 1] 2 - Aa 2 M(d - Z)^- 1 
[(d - l)Mu d - 3 - 1] x [a 2 u - (d - 3)M u d ~ 4 ] 



-2a 2 u [(5 - d)M- 



.d-3 



1] 



+ [a 2 u -(d- 3)Mu 



d-41 



(B7) 



The discriminant associated to this equation is given by 
16Ma 2 (d — 3)u d_1 A 2 , where we have introduced A„ = 
a 2 u 2 — 2Mu d ~ 3 + 1. In order to write the solutions of 
equation (|B7|) , it is convenient to consider the expression 

[(d-l)Mu d - 3 - l] 2 - Aa 2 M(d-Z)u d - 1 = Q_Q H 

where 



(B8) 



Q± = l-(d- l)Mu d - 3 ± 2a^J(d - 2,)Mu d - x . (B9) 

Then the solutions of Eq. (|B7j) can be written as 

1 



2 2 Q±A U - Q^ 

x u = 



Q+Q- Q T 

As an alternative, we can use the identity 



(A„ - Q T ) . (BIO) 



A„ - Q T = 



aV"± \ (d-3)Mu d ~ 



to cast the solution for x in the simple form 



<V^± y/(d- 3)Mu d ~ 



(Bll) 



(B12) 



where the upper sign in the foregoing equations applies 
to counterrotating orbits, while the lower sign applies to 
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corotating orbits. Inserting expression (|B12|) in equation 
(|B5I) and using the relation L = aE + x, we get the 
following expression for E, 



E = 



1 - 2Mu d ~ A T aJ{d - 3)M 



,d-l 



(B13) 



and the angular momentum associated with the circular 
geodesies, 



L = T 



yj(d - 3)Mu d - 
\/uQ T 



1 + a 2 u 2 ± 2a\ 



d- 3 



(B14) 

In order to investigate the stability of circular timelike 
orbits we must compute the second derivative of V with 
respect to r for the values of E and L specific to circular 
orbits. Differentiating Eq. (|B1[) we find 



V" = \2r z (E A -l) + 2{d-4)(d-5)Mr 3 - d x 2 



-4aEx - 2a 2 + 2(d - 6)(d - 7)Mr 



5-d 



- 2x z 
(B15) 



Substituting x and E from Eqs. (|BT2|) and (|Bl"3|) , the 

above expression for V" becomes 



V" 



2(d-3)Mu 



d-5 



2{d-l)Mu d - 6 + (d-5) 



±8ayJ(d~3)Mu d - 1 + (d- l)a 2 u 2 



(B16) 



The term within square brackets in the foregoing equa- 
tion is equal to (d — 1)A« — 4Q T , so that the second 
derivative of V reduces to 



4Q= 



(B17) 



This expression shows an explicit dependence on the 
spacetime dimensionality d. To analyze the sign of V" 
in Eq. (|B17p it will be helpful to distinguish between dif- 
ferent values of d. Since E, L and x = L — aE must be 
real, the functions A u and Q± are such that 



A u > Q± > 0. 
For d > 5 the above conditions lead to 



(d-l)A u >4Q T 



V" > . 



(B18) 



(B19) 



This means that there are no stable timelike circular or- 
bits for spacetimes with d > 5. This generalizes previous 
work by Tangherlini on non-rotating higher-dimensional 
black holes [2y] and by Frolov and Stojkovic on five- 
dimensional rotating black holes [4^ . 



2. The orbital frequency of circular geodesies 

The orbital frequency il = dip/ 1 dt associated to circular 
timelike geodesies is given by 



(L - 2Mu d ~ 3 x) u 2 
(1 + a 2 u 2 ) E - 2aMu d ~ 1 x 



(B20) 



The foregoing expression can be simplified by considering 
the following identities: 



, , J(d - 3)Mu d - 4 . 
L — 2Mu x — =F — [ — A„, (B21) 

y / uQ T 

(l + a 2 u 2 ) E - 2aMu d ' 1 x 

A " ' _ " ' (B22) 



l=Fa\/(d-3)Mu d - 1 



Substituting (|B2l|) and (|B22f into equation ([B20)> . we 
obtain 



= 



T v / (<i-3)Afu' i - 



lTa- v /(^-3)Mu d - 1 



(B23) 



By considering equation ([79| for D c and the relation 
il c = 1/-D C , we find a similar expression for circular null 
geodesies: 



W(d-3)M- 



,d-l 



(B24) 



l Ta y'(d-3)Afu^ 1 
Consequently, the ratio fi/f2 c varies with r as follows: 



(B25) 



APPENDIX C: SCHWARZSCHILD-ANTI-DE 
SITTER SPACETIMES 

In this Appendix we consider Schwarzschild-anti-dc 
Sitter spacetimes, and we show that the analogy between 
unstable circular orbits and black hole QNMs is not triv- 
ially extended to non-asymptotically flat backgrounds. A 
higher-dimensional Schwarzschild-anti-de Sitter solution 
is a solution of 



(rf-2)(rf-l) 

^ab H ^7 9ab — U , 



2L 



ads 



(Cl) 



with cosmological constant A = — (d — 2)(d — l)/(2i a d s ) 
and typical curvature radius L a ds- We consider 
the simplest black hole solution: the d-dimensional 
Schwarzschild-anti-de Sitter solution. The line element 
is given by Eq. (f22|) with 



f(r) = g(r) = 



T 2 



1 - 



„d— 3 



„d-3 



(C2) 



The quantity r is related to the mass M of the space- 
time, 



M 



(d - 2)A d - 2 r d 
16tt 



d-3 



(C3) 



and the horizon radius r + is the largest real root of f(r) 
0. 
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For circular null geodesies, an exact solution can be 
found with 



2^ 



(d-1) 1 + 



T2 

ads 



r+ ■ 



(C4) 



Since V"(r c ) = L 2 (2d — 6) r c 4 , circular null geodesies 
are unstable. The angular velocity at this radius is given 
by 



(d-3) 



1 



d - l 
d-3 



1 



r 2 



ads 



(C5) 

and reduces to Cl c » 1/L a d s for large (r + /L a & s 3> 1) 
black holes. The calculation of the instability exponents 
proceeds trivially. As in the case of higher dimensional 
Schwarzschild black holes, Eq. ([50]) . we find 



x/n c = Vd- 3. 



(C6) 



If our main result, Eq. (|46|) . were valid in asymptoti- 
cally AdS spacetimes, in the eikonal limit we would get 
wqnm/^c = l — iy/d — 3 (n+1/2). However, according to 
both analytical 56] and numerical results [59, 60, 61, 62], 
the imaginary part of QNM frequencies in this back- 
ground increases monotonically with increasing r + . This 
dependence on r + cannot be explained by the circular 
null geodesic analogy. In hindsight, this failure is not 
too surprising. AdS spacetimes are not globally hyper- 
bolic; boundary conditions at infinity must be taken into 
account. On the other hand geodesic calculations are 
local, and they carry no information about spatial in- 
finity. The available analytical and numerical results 
[Hrj [59l loOl . f6lL [621 ] indicate that the damping timescale 
is smaller than indicated by the geodesic calculation, 
Eq. (|6"2")) . Perhaps the disagreement could be explained 
by arguing that null particles reach spatial infinity on 
timescales faster than the geodesic timescale, and there- 
fore one would have to correct for this. 
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